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ABSTRACT 

We compute the evolution of the space-dependent mass distribution of galaxies 
in clusters due to binary aggregations by solving a space-dependent Smoluchowski 
equation. From the solutions we derive the distribution of intergalactic distance for 
different ranges of mass (and of corresponding magnitude). We compare the results 
with the observed distributions, and find that the different degrees of luminosity 
segregation observed in clusters are well accounted for by our merging model. In 
addition, the presence of luminosity segregation is related to dynamical effects which 
also show up in different but connected observables, such galaxy velocity profiles 
decreasing toward the center and X-ray measured /3-parameters smaller than 1. We 
predict both luminosity segregation and the observables above (being a product of 
binary aggregations) to be inversely correlated with the core radius and with the 
galaxy velocity dispersion; we discuss how the whole set of predictions compares with 
up-to-date observations. 

Subject headings: galaxies: clustering - galaxies: intergalactic medium - galaxies: 
X-rays - hydrodynamics 



1. Introduction 

The dynamical evolution of galaxy clusters is currently believed to go through two major 
phases: in the first, usually referred to as violent relaxation (Lynden-Bell 1967), the evolution 
is controlled by a collective potential and results in a Maxwell velocity distribution of galaxies; 
in the second, the dynamics is dominated by two-body processes, and binary (both elastic and 
inelastic) collisions drive the evolution. In fact, in this latter phase, for ordinary galaxy sizes and 
separations the collision time scale is much less than the Hubble time. 

Although a complete theoretical description of the two-body phase of dynamical evolution 
of clusters is still lacking, observations, N-body simulations and computation based on statistical 
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methods (Monte Carlo and Fokker-Plank simulations) have concurred in enlightening many 
dynamical properties of clusters in this stage. Such properties are characterized by a large cluster- 
to-cluster variance, and include the following: the presence of a velocity bias b 2 , = (v 2 )/o~ 2 < 1 of 
the galaxy velocity dispersion (v 2 ) 1 / 2 with respect to the dark matter's a (see N-body simulations 
by Carlberg k Dubinski 1991; Evrard, Summers k David 1994; Katz k White 1993; Carlberg 
1994; Summers, Davis k Evrard 1995); galaxy velocity dispersion profiles decreasing toward the 
cluster center (see observations by Kent k Sargent 1983; Sharpies, Ellis, k Gray 1988; Girardi 
et al. 1996); mass segregation, i.e., the tendency of more massive galaxies to be located near the 
cluster center (see simulations by Roos k Aarseth 1982; Farouki, Hoffman k Spencer 1983) with 
the associated luminosity segregation (observed in several clusters, see Rood et al. 1972; Oemler 
1974; White 1977; Dressier 1978; Quintana 1979; Sarazin 1980; Kent k Gunn 1982; Oegerle, 
Hoessel k Ernst 1986; Binggeli, Tammann k Sandage 1987; Dominguez Tenreiro k Del Pozo-Sanz 
1992; Stein 1996). 

Due to the large variance observed in the above effects, appreciable uncertainties exist about 
their dependence on the characteristics of the system, although a general trend in the sense of 
larger effects for smaller clusters might be inferred from the data, when only relaxed cluster (with 
no prominent substructure or asphericity) are considered. 

A complete description of the post-virialization phase of galaxy cluster should be able to 
connect all the above effects and to explain the observed variance in terms of dynamical properties 
of clusters. In previous papers (Fusco-Femiano k Menci 1995, hereafter Paper I; Menci k 
Fusco-Femiano 1996, hereafter Paper II) we showed that the loss of kinetic energy in inelastic 
galaxy collisions (binary aggregations) in clusters with a ^ 900 km/s can significantly change the 
velocity distribution of galaxies. The model not only simultaneously accounts for the velocity bias 
and for centrally decreasing velocity dispersion profiles, but also predicts a correlation of such 
effects with the shape (the core radius) and the depth (the dark matter velocity dispersion) of the 
cluster potential wells in good agreement with observations. In addition, the aggregation model 
can succesfully connect the above dynamical effects to other observed properties of galaxy clusters 
(Paper I), such as the /3-parameter (expressing the ratio of the galaxy orbital specific energy to 
the specific energy of the X-ray emitting plasma) and the Butcher-Oemler effect (see Cavaliere k 
Menci 1993). 

Thus, aggregations seem to constitute a leading mechanism in the post-virialization phase of 
clusters with velocity dispersion ^ 900 km/s, while in larger clusters, they are highly suppressed 
due to the large galaxy relative velocities (see numerical results by Richstone k Malmuth 1993). 
To further assess the role of galaxy merging in the two-body dynamical phase, we address here 
the problem of luminosity segregation (hereafter LS). This is expected to be generated when 
aggregations are effective, since merging builds up larger galaxies mainly in the central regions, 
where the larger density favours binary aggregation. Thus, we extend our treatement of galaxy 
inelastic collisions, based on the solution of a collisional Boltzmann-Liuoville equation, to include 
position-dependent mass spectra of interacting galaxies. The prediction of our model will be 
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compared with observational results, focussing on the correlation of the segregation effects with the 
properties of the clusters such as the richness, the density distribution and the velocity dispersion. 
Finally, we shall show how the merging model connects LS with the dynamical effects discussed 
above and with their X-ray counterparts. 

The paper is organized as follows. In sect. 2 we discuss the collisional Boltzmann equation 
and describe our solutions for the evolution of the position-dependent mass distribution. In 
Sect. 3a we decribe the standard method we use for the comparison with data, based on the 
auto-correlation function for galaxies of different mass. The comparison is performed in Sect. 3b. 
Sect. 4 is devoted to discussion and conclusion. 



2. Radius-dependent mass distribution from binary aggregations 

2.1. The Boltzman Equation for Merging Galaxies 

The evolution with time t of the distribution f t (M,r,v) of interacting galaxies with velocity 
v and mass M at the position r inside the gravitational potential ip of a cluster can be described 
by the collisional Boltzmann equation. Assuming spherical symmetry, the latter can be written in 
spherical coordinates r = (r,9,<f>) and v = (v r ,V0,v,p) for the distribution f(M,r,v r ,Vt) as follows 
(see, e.g., Saslaw 1985) 



o-zw \ i df t (M,r,v r ,v t ) /vf dip\ v r vf df t (M,r,v r ,v t ) 
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where v\ = v 2 + v 2 ^ is the square tangential component of velocity, and the gravitational cross 
section for interactions E depends on the relative velocity v re ; = v' — v". The velocity v" in the 
first integral in eq. (1) is related to v and to v' by the requirement of momentum conservation 
M' v' + (M — M') v" = M v. Here we have assumed that the galaxies do not gain or loose mass 
via processes other than merging. To obtain a fully self-consistent description, eq. (1) should be 
complemented with the Poisson equation for the gravitational potential ip with a source term 
/ dM dv r dv 2 f t (M,r,v r ,v t ), which includes the galaxy distribution itself. However, since we are 
interested in the post-virialization phase of cluster evolution where the potential is essentially fixed, 
we shall assume a King potential i/j(r) and follow the evolution of the galaxy mass distribution at 
different radii. 

Since our aim is to probe the effectiveness of interaction in producing mass segregation, we 
indroduce some approximations (a discussion on them is given in the final section). First, we 
assume the velocity distribution to be independent on the mass and on the spatial distribution of 
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galaxies, so that the distribution in eq. (1) can be factorized into a velocity distribution p(v) and 
a position-dependent mass distribution N(M,r,t). Such approximation does not actually hold 
(see Paper I) but, as we discuss in the final Section, for our purpose in the present paper it is a 
conservative assumption. Second, we assume that the radial and tangential velocity distribution 
are mutually independent and both normally distributed. In this case, integration of eq. (1) over 
velocities leads to the following position-dependent Smoluchowski equation 

d t N(M,r,t) = -JdM'N(M',r,t)N(M-M',r,t)(£(M',M-M')vrei) + 

dM' N(M, r, t) N(M', r, t) (E(M, M 1 ) v rel ) (2) 
where the average () is over the velocity distribution. 

r 2-i 

The cross section is given by (Saslaw 1985): S(M, M') = e(v re i/v g ) it (r 2 + r' 2 ) 1 + where 

r and r' are the radii of the interacting galaxies (proportional to M 2 / 3 ) and v g oc G(M + M')/R is 
the escape velocity at closest approach R ~ (r + r'). The efficiency e is determined from N-body 
results (see Richstone &; Malmuth 1993) and is zero when v re i ^ 3v g , so that aggregations are 
highly suppressed in very rich clusters. It is convenient to express all quantities in terms of the 
adimensional mass m = M/M* normalized to the characteristic mass M* (that corresponds to a 
galaxy with characteristic luminosity L*). From r ~ (M/p) 1 / 3 , the relation r 2 = r 2 ^ m 2 / 3 follows. 
Then, the cross section reads 

S(m, m') = e(t; rc ,/ Vff ) Trr 2 , (m 2 / 3 + m' 2 / 3 ) [l + (m 2 / 3 + m' 2 / 3 ) u^/^] (3) 

where r s * and "U^r^ axe the radius and the 3-D internal velocity dispersion of a galaxy, 
respectively. 



2.2. Initial Conditions 

We assume the galaxy distribution to be initially (after the cluster formation and virialization) 
factorized in a mass distribution P(m) times a King spatial profile, which will be subsequently 
mixed up by the 2-body dynamical evolution. Then 

MM,r^o= (1+ ra ; 2)3/2 P(m), (4) 

where we take for P(m) the Press & Shechter shape P(m) = m a ~ 2 e _fe2 ^ m2a I 2 (the index a 
depends on the spectrum of cosmological perturbations and is in the range - 0.3 at the scale of 
galaxy clusters) and x = r/r c is the distance from the cluster center in units of the core radius r c 
of the King profile. The constant n Q is taken as to yield the total number N to t of galaxies inside 
the cluster virial radius R v (from the virial theorem R v = GM/3a); thus 

120 = ~A 3 T T ' ( 5 

4 vr r% I R I m 
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where Ir = Jo*" dxx 2 /(l + x 2 ) 3 / 2 and 1m = Jo°° dmP(m) are the adimensional integrals of the 
initial spatial and mass distributions, respectively. 



2.3. Numerical Solutions 

To integrate eq. (2) we first write it in a completely adimensional form for the normalized 
r-dependent mass distribution rit(m,r) = N(m,r,t)/n . We define the adimensional time 
variable r = t/t cr ~ 2 10 9 yr [R v /1 Mpc] [cr/10 3 km/s] -1 in terms of the the cluster crossing time 
t cr = 2R v /ct. The adimensional velocities v = v/a are normalized to the dark matter velocity 
dispersion a. The corresponding adimensional interaction rate r](m, to') = n Q £ R v v can be 
computed from eq. (3) and (4). Then, the Smoluchowski eq. for the normalized mass distribution 
n T (m, r) = N(M,r,T)/n can be recast in the form 

1 f m 

d T n T (m,r) = - dm n T {m' ,r) n T {m — m! ,r) (rj(m' ,m — m 1 )) + 



2 



o 



dm'n T (m, r) n T (m', r) (r](m, to')) (6a) 



V (m, to') = I % Pf- v rel (m 2 / 3 + to' 2 / 3 ) [l + (to 2 / 3 + m' 2 / 3 ) v 2 g Jv 2 rel ] . (6b) 

From the form (6) it is evident how (for constant R v /r c ratio) the effect of aggregation is larger 
for clusters with small core radius (galaxies in the center are more concentrated) and with a larger 
number of galaxies Ntot- 

The average of the aggregation rate in eq. (6b) 

r r\vi-v 2 \=3v g 

(rj) = da dviv 1 p(v 1 )dv2V 2 p(v 2 )'n{v 1 - v 2 \/v g ) , (7) 



is over the velocities v\ and v 2 (normalized to the dark matter velocity dispersion a) of galaxies 
colliding with relative angle a; the condition |vi — V2I < 3v g accounts for the efficiency e(v re i/v g ). 
We assume the distribution of velocities p(v) = (l/27r) -3 / 2 e~ v / 2 to be Gaussian, as expected after 
violent relaxation (Lynden-Bell 1967). Note that, for clusters with a <^ 900 km/s, eq. (7) yields 
significant averaged aggregation rates (rj) assuming a 3-D internal velocity dispersion v g * = 300 
km/s for an L* galaxy with r gif = 60/i _1 kpc Q. The adopted value of v g * correspond to a circular 
velocity of ~ 220 km/s; such value is consistent with that derived from the Faber- Jackson relation 
for an L* galaxy, and with the measurements by Tonry & Davis (1981); Dressier (1984); Dressier 
(1987). The adopted value of r gif (which refers to the dark halo of an L* galaxy) is a conservative 
one, when compared with observational results from absorption lines, measured by Steidel (1995); 
Lanzetta et al. (1995); Barcons, Lanzetta & Webb (1995). 



1 In the text we adopt h = 0.5 for the Hubble constant H — 100 h km/s/Mpc. 
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Equation (6) is integrated up to r = 5 with time increments At = 1/500 and mass step 
Am = 1/500, from a minimum mass m,i n f = 10~ 2 to a maximum mass m sup = 10 2 (integrating 
up to larger times do not affect sensitively our results). When aggregations are effective, the final 
mass distribution will be changed from the initial one only in the central core where the galaxy 
density is larger and binary aggregations are favoured. Thus, in the core larger galaxies will form 
via binary merging, while the initial mass distribution remains unchanged in the outer regions. 
The evolution of the mass distribution at different radii in a typical cluster (with N to t = 1000, 
a = 800 km/s and r c = 250 h^ 1 kpc) is shown in figure 1. The distribution flattens in the central 
region due to the disappearence of small galaxies which aggregate to form larger ones. Since 
aggregations between galaxies cause a loss of orbital kinetic energy, we expect such effect to be 
correlated with galaxy velocity dispersions smaller in the central regions, i.e., with velocity profiles 
falling toward the center (as we discussed in Paper II); this is actually the case, as is shown 
in figure 2. A further effect is that brighter galaxies (which form from mainly in the central, 
denser regions) will have smaller relative separations. This latter effect is that observed in several 
clusters, as we discuss in the next session. 

The strenght of the above effects depends on the cluster properties, which, in our model, enter 
only through Ntot, cr and r c as is shown by eqs. (6). E.g., for given Not and <r, in clusters with 
large r c merging will be less effective (see the merging rate in eq. 6b) because the total number of 
galaxies is spread out in a larger region. 

3. Comparison with Observations 
3.1. Method 

In the literature (Capelato et al. 1980; Dominguez-Tenreiro & del Pozo-Sanz 1988) the LS 
has been quantified in terms of the cross correlation function 

U a (s)= f #d 2 rn a (r)n a (r + s) (8) 
Jv 

between densities of galaxies separated by a distance s in a given magnitude range [a] ((f) is the 
angle between r and s), in a region V. The distance distribution function for pairs in a given class 
is then given by 

P a (s)ds = 2vrsdsn a (s) . (9) 

If the position of the peak in the distribution P(s) changes depending on the magnitude class [a] 
a LS is present. The average separation of galaxies in the magnitude class [a] derived from eq. (8) 
reads 

P Sraax P Sraax 

A„= / dssP a (s)/ / dsP a ( S ) , (10) 

JO JO 

where S max is the maximum intergalactic distance. Galaxies belonging to a class [a] will be called 
"segregated" with respect to those in the class [a'] if A a < \ a > . 
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Such an effect has been measured in several clusters (Capelato et al. 1980; Dominguez- 
Tenreiro & del Pozo-Sanz 1988; Yepes, Dominguez-Tenreiro & del Pozo-Sanz 1991; Yepes & 
Dominguez-Tenreiro 1992) and a fit to the observed P(s) for different classes of magnitudes has 
been given by the same authors. We consider all the clusters where such an analysis has been 
performed, except for Perseus, whose prominent substructure (Gallagher, Han & Wyse 1996)) 
makes it too complex to be describable in the framework of our model. They are reported in 
Table 1 together with their core radius, velocity dispersion and observed number of galaxies inside 
a distance R ma x- In the same table are listed the magnitude ranges [a] for which we compute the 
correlation function. 

To compare with such observational material we proceed as follows: 

• For each observed cluster, we compute the number of galaxies N to t enclosed inside Ry 
for the whole mass range 0.01 < M/M* < 10 (used in the numerical integration of eq. 6), 
corresponding to a luminosity range 0.01 < (L/L*) 7 < 10 for a M/L oc L 7_1 (here we take 7 = 1, 
but see discussion in Sect. 4 for the effect of changing M/L). 

In practice, Ntot is computed extrapolating the observed number N Q bs of galaxies (inside a 
radius R ma xi see Table 1) both in space (up to R v , using a King profile) and in luminosity (for the 
whole luminosity range discussed above, using a Shechter luminosity function). 

The resulting N tot is given as an input for the solution of eq. (6). together with the cluster 
core radius r c and the dark matter velocity dispersion a. The latter is derived from observed 
galaxy velocity dispersions (assuming no velocity bias) or, when the latter are not available, from 
the X-ray temperature T = (finiH /k j3)a 2 (when measures of (3 are not available, we shall assume 
(3 = 1). The resulting values (with the references to the corresponding observations) are given 
in Table 1. The reported a are affected by uncertainties Aa/a < 20% due to intrinsic errors in 
the measurements of velocity dispersions or X-ray temperatures and (when the estimate of a is 
obtained from T with no available measurements of (3) to the indetermination of (3. However, we 
stress that errors in a (as well as those on N b s ) do not affect sensitively the LS effects resulting 
from our model. A quantitative discussion of the effect of variations of all the input parameters is 
given in Sect. 3.3. 

• For each cluster the r-dependent mass distribution is found integrating numerically eq. (6). 

• We divide the computed mass distribution at each radius according to the classes of 
apparent magnitudes (see Table 1) which have been used in the analysis by Capelato et al. 1980; 
Dominguez-Tenreiro & del Pozo-Sanz 1988; Yepes et al. 1991; Yepes & Dominguez-Tenreiro 1992. 
To pass from magnitude to mass ranges we use the M/L ratio discussed above. 

• We compute the distance distribution P(s) resulting from our model and compare it with 
the fit to observational results found in the literature. For each cluster, the average separation 
X a corresponding to P a (s) is computed for all the magnitude classes [a] and compared with the 
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observed values. 

When the cluster characteristics are such as to make aggregations effective, larger galaxies 
form preferentially in the central, denser regions (see Sect. 2 and fig. 1) where the intergalactic 
separation are smaller. In this case, the distributions P a (s) will be peaked at smaller separations 
for brighter magnitude ranges [a]. Such shift of the peak with [a] can be expressed by the ratio 

bx = A1/A3 (11) 

of the average distances of the brightest class to that of the faintest class. 



3.2. Results 

The distributions P a (s) for the different clusters are shown in fig. 3. The ratio b\ is shown 
in Table 2 for our predictions and for the corresponding observations. The agreement with 
observations is remarkable. The different degrees of segregation (expressed by values b\ < 1) 
observed in the sample is well accounted for by our merging model, and is directly related to the 
cluster characteristics as follows: for a given total number of galaxies N to t, clusters with small 
core radii have denser central regions, so that aggregations are more effective and segregation 
is enhanced. The lack of LS in A2111 is explained in terms of large core radius coupled with a 
limited total number of galaxies. As a confirm to such an interpretation, we observe that the more 
pronounced segregation takes place in A2670 which is characterized by the smallest core radius in 
the sample. However, LS can occur also in clusters with large r c if the total number of galaxies 
Nt t is large enough or if a is very low. In fact, the cluster 0004.8-3450 shows a significant LS due 
to Ntot — 2000, while the segregation in the Fornax cluster is mostly due to its very low velocity 
dispersion cr=320 km/s. 

Note the peculiarity of cluster A2218. The model is in very good agreement for what 
concerns the observed LS of the two brightest magnitude classes with respect to the third one. 
However, the real data show that the two brightest classes are characterized by an anti-segregation 
between them, which is not accounted for by our model. We attribute such mismatch to 
substructures/anisotropies which our model (based on the schematic assumption of isotropy) 
cannot reproduce. In fact, recent analysis (Squires et al. 1996) of A2218 indicates the presence of 
collision of subclumps, with associated elongated structures in the plasma disposition. 

As observed above, our model predicts the aggregating galaxies to loose part of their kinetic 
energy. The brightest galaxies in a cluster showing segregation are then expected to have velocity 
dispersion profiles decreasing toward the center, where the larger density favours aggregations. 
The computed result for the galaxies belonging to the magnitude classes 1 and 2 in A2670 (see 
figure 4) confirm such expectation and are consistent with the available data for such cluster 
(Sharpies, Ellis, & Gray 1988; see also Yepes & Dominguez-Tenreiro 1992). The same calculation 
for A2111 (see figure 4) shows no positive gradient in the profiles, which is the counterpart of the 
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lack of luminosity segregation. Actually, both effects are tightly connected in our model. 

3.3. Varying the Input Parameters 

Here we discuss the effects of variations of the input parameters with respect to the reference 
values in Table 1. The segregation parameter b\ decreases (indicating larger segregation) for 
increasing N b s , and for decreasing a and r c (i.e., for increasing merging efficiency). However the 
variations with N b s and a are very mild. This makes our results robust with respect to the errors 
associated to those parameters: a 20 % error in AT o6s or in a results in Ab\/b\ < 3%. The errors 
in r c are more important: Ar c /r c = 20% yields Ab\/b\ < 12%. 

The results for LS do depend on the M/L ratio which for the sake of simplicity, we assumed 
to be constant. However, the main results presented here holds also for M/L oc L 7_1 with 
3/4 < 7 < 4/3. This is illustrated in figure 5, where we show the distance distribution functions 
(for the parameters of cluster A2670) derived from the same dynamics (i.e., with the same mass 
segregation) but with 7 = 3/4 and 7 = 4/3. 

4. Conclusions 

We have shown that a detailed model for the dynamics of galaxies aggregating in the potential 
wells of clusters predicts luminosity segregation (LS) effects of the kind observed in real clusters. 
The correlation of the strenght of the effect with the properties of the clusters predicted in our 
model is in agreement (see figure 3 and Table 2) with that observed for the limited sample of 
clusters for which LS has been subject to accurate quantitative measurements. In particular, we 
predict the effectiveness of aggregations, and hence the degree of LS, to be directly correlated with 
the number of galaxies in the cluster and inversely correlated with the core radius and with the 
velocity dispersion (see eq. 6b). 

The results do not depend on the detail of the initial mass distribution of galaxies in 
clusters, which we assume to have a Press & Shechter form with spectral parameter a = —2; such 
independency is due to the properties of the asymptotic solution of the Smoluchowski equation 
(describing the evolution of the position-dependent galaxy mass function in our model) and can 
be traced back to the non-linear nature of such equation. 

Our results are robust with respect to uncertainties in the input quantities and to the adopted 
L(M) , as shown in Sect. 3.3. As for our assumption of fixed galaxy velocity distribution, this 
does not hold when aggregations are effective (see Paper I and II). However, the merging-induced 
shift of ~ (10 — 15)% of the velocity dispersion toward smaller values increases the efficiency of 
aggregations, so that our assumption is actually conservative. We have re-run our computation 
for shifted velocity distributions and found results almost indistinguishable from those presented 
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here. Finally, we stress that no attempt of parameter optimization has been performed. An even 
better agreement could be found if the cluster parameters were suitable tuned. 

As for the big picture of the evolution of cluster in the two-body dynamical phase, our model 
focus on the effects of inelastic collisions not considered in previous works on this subject. In 
particular, the Fokker-Plank approach by Yepes & Dominguez-Tenreiro (1992) considers only 
elastic collisions by means of a "mean field" approximation with the input parameters chosen 
from a grid of models to show that, within the set of models, it is possible to match the observed 
segregation effects. 

Here we solve the collisional Boltzmann equation including inelastic collisions, using, for the 
input model parameters, the measured values. Though the latter are subject to errors, we showed 
(in Sect. 3.3) that the model is robust to uncertainties < 12% in the parameters. Our results show 
that inelatic collisions produce appreciable dynamical effects for clusters with one-dimensional 
velocity dispersion ^ 900 km/s. Such effects show up in different but connected observables: the 
velocity bias (due to the average loss of kinetic energy in inelastic collisions) b v « 0.8 — 0.9 can be 
observed in X-rays in the form of /3-parameter (see Cavaliere & Fusco-Femiano 1976) (3 = b 2 v < 1; 
centrally rising velocity dispersion profiles (due to the differential loss of kinetic energy at different 
radii) are now being measured with great accuracy in the optical (Girardi et al. 1996); different 
average separations of massive galaxies with respect to the faint ones b\ ~ 0.8 — 0.9 (see eq. 5), 
i.e., luminosity segregation (due to the differential mass growth from aggregations at different 
radii) have been measured in different clusters (see references cited in this paper). 

We note that, when interpreted in terms of merging-driven evolution, all the above effects are 
predicted to have the same dependence on the cluster parameters, i.e., to be larger for clusters with 
smaller core radii r c and galaxy velocity dispersions a, although the strength of the c-dependence 
is mild for LS effects. 

The observational tests for such predictions are critically affected by the presence of clusters 
with anisotropics and/or substructures in the observational sample. An inverse correlations of the 
(3 < 1-effect with r c has been found by, e.g., Jones & Forman 1984, while the anti-correlation with 
a has been pointed out by Kriss et al. (1983); Jones &: Forman (1984); Edge & Stewart (1991); 
Bird, Mushotzky k, Metzler (1995); Jones et al. (1997) but has not been confirmed by the analysis 
by Lubin & Bahcall (1993) and by Girardi et al. (1996). 

For velocity dispersion profiles decreasing toward the center, the observational situation is still 
unclear. An anti-correlation with a has been inferred (see Paper II) from the analysis by Girardi 
et al. (1996) of a sample of 37 clusters, when clusters with prominent substructures are excluded; 
however, the detection of such correlation from the data (see also den Hartog, & Katgert, 1996) is 
made difficult by the presence of anisotropies and/or substructures, which can hurt or destroy the 
effect of the inelastic collision in the two-body relaxation phase. 

As for the LS, the mild (inverse) dependence on a of the LS effect from merging makes 
difficult to observe such correlation. However, in our model, the strong inverse correlation of LS 
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with r c predicted by our model is confirmed the data analysis by Yepes, Dominguez-Tenreiro & 
Del Pozo-Sanz (1991) on the very limited sample of clusters. LS data for a larger sample of cluster 
with measured r c would definitely clarify the issue. 

Finally, we note that the correlations between different but connected observables predicted by 
the aggregation model makes it testable already at the present stage of observational capabilities. 
Further observational progress (in particular in measuring in detail velocity dispersion profiles and 
X-ray temperatures) can definitely probe the predictions of the merging picture, thus assessing 
the role of aggregations in the dynamical evolution of clusters. 

We thank the referee for keen suggestions and comments. 
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Fig. 1. — The position dependent mass distribution obtained from the solution of eqs. 6, for a 
cluster with the parameters given in the text (arbitrary units are used for r and M). Top panel 
refers to the initial condition, while the bottom panel to the evolved distribution at t = 5t cr . Note 
the flattening at the low-mass end in the central region 

Fig. 2. — The corresponding integrated velocity dispersion calculated as shown in Paper II. Note 
the decrease toward the central region, corresponding to the loss of orbital energy by aggregations 
occurring in the cluster center 

Fig. 3. — The computed intergalactic distance distribution for the clusters with parameters listed 
in Table 2. The solid curve refers to the brighter magnitude class, while the dotted line to the 
fainter. All the curves have been computed for a constant M/L 

Fig. 4. — The integrated velocity dispersion profiles computed as shown in Paper II, for the cluster 
A2670 and A2111 

Fig. 5. — The computed intergalactic distance distribution for the cluster A2670, for different values 
of the M/L oc L 7_1 (see text) ratio. Top panel refers to 7 = 3/4, while bottom panel to 7 = 4/3. 



TABLE 1 
Analyzed Clusters of Galaxies 



Cluster Name 


z 


r c 


Fornax 


N ois a 




Group 


Magnitude Range 


A1758 


0.28 


0.72 


1.25 


320 


1000 


1 


18.09 < m F <20.10 


A2111 


0.23 


0.94 


1.08 


269 


970 


2 
3 
1 


20.10< m F <20.90 
20.90< m F <21.70 
17.43< m F <19.43 


A2218 


0.171 


0.43 


0.87 


306 


800 


z 

3 
1 


10 /iQ^ ™ /on 01 
iy.4o<. Tflp <.2V.21 

20.21< m F <21.00 
17.12< m ; <19.12 


A2670 


0.076 


0.198 


0.59 


220 


600 


2 
3 
1 


19.12< m f <20.15 
20.15< m s <21.18 
15.40< m 6 <18.00 


Coma 


0.023 


0.3 


1.3 


400 


810 


2 
3 
1 


18.00< m b <19.00 
19.00< ra s <20.00 
11.80< m 25 <14.50 


Fornax 


0.005 


0.36 


0.73 


68 


320 


2 
3 
4 
1 


14.50< m 25 <15.50 
15.50< m 25 <16.50 
16.50< m 25 <17.50 
10.20< m B <13.50 


0004.8-3450 


0.114 


1.20 


1.87 


333 


805 


2 
3 
1 


13.50< m B <15.50 
15.50< m B <16.50 
16.32< m B <19.00 



2 19.00< m B <20.00 

3 20.00< m B <21.00 



a Total number of plate galaxies in the circular region of radius R m ax 

b The values in the Table, used as inputs in our model, are derived as described in the text from the following 
observational material. 

A1758: T = 7 keV and T = 8 keV (two clumps) from Ebeling et al. (1996), corresponding to a = 1000 — 1200 km/s 
assuming (3=1. 

A2111: T = 6 KeV from Wang k Stocke (1993) yields a = 970 km/s assuming /? = 1. 

A2218: direct estimates from Girardi et al. (1997) give a = 700 - 800 km/s. 

A2670: T = 3.9 keV from David et al. (1993) yield a = 600 for j3 = 0.6 (Jones & Forman 1984). 

Coma, Fornax and 004.8-3450 (also refferred to as A2721): direct estimates from Fadda et al. (1996). 

c Original references for photometric data are: Butcher, Oemler & Wells (1983) for A1758, A2111, A2218; Sharpies, 
Ellis, k Grey (1988) for A2670; Ferguson (1988) for Fornax; Carter (1980) for 0004.8-3450; Godwin k Peach (1977) for 
Coma. 
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TABLE 2 
Lenght Scales 



Cluster Name b\ (observed) b\ (model) 



A1758 


0.80 


0.78 


A2111 


0.94 


0.96 


A2218 


0.85 


0.78 


A2670 


0.55 


0.56 


Coma 


0.83 


0.87 


Fornax 


0.88 


0.89 


0004.8-3450 


0.85 


0.77 



Errors in the observed lenght scales are < 20%. For the errors in 
the the model values, see Sect. 3.3. 
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